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ABSTRACT 

We present new H and K bands spectroscopy of 15 high luminosity active galactic nuclei (AGNs) at red- 
shifts 2.3-3.4 obtained on Gemini South. We combined the data with spectra of additional 29 high-luminosity 
sources to obtain a sample with 10 452 < XLx (5100A) < 10 47,3 ergs s _1 and black hole (BH) mass range, using 
reverberation mapping relationships based on the H/3 method, of 10 88 — 10 10 7 M @ . We do not find a corre- 
lation of L/Leaa with Mbh but find a correlation with AL^(5100A) which might be due to selection effects. 
The L/Ledd distribution is broad and covers the range ~0.07-l .6, similar to what is observed in lower redshift, 
lower luminosity AGNs. We suggest that this consistently measured and calibrated sample gives the best repre- 
sentation of L/Leaa at those redshifts and note potential discrepancies with recent theoretical and observational 
studies. The lower accretion rates are not in accord with growth scenarios for BHs at such redshifts and the 
growth times of many of the sources are longer than the age of the universe at the corresponding epochs. This 
suggests earlier episodes of faster growth at z >, 3 for those sources. The use of the C IV X 1549 method gives 
considerably different results and a larger scatter; this method seems to be a poor Mbh and L/Lgdd estimator at 
very high luminosity. 

Subject headings: Galaxies: Active - Galaxies: Nuclei - Galaxies: Quasars: Emission Lines 



1. INTRODUCTION 

Studies of black-hole (BH) growth at various redshifts, and 
the comparison with galaxy evolution and star formation, has 
been a very active area of research for several years. In partic- 
ular, there are several suggestions that very massive BHs grew 
faster at early epochs while the growth of less massive BHs 
extends over longer periods and is significant even at z — 0. 
For example, Marconi et al. (2004) used the X-ray luminosity 
function of active galactic nuclei (AGNs) to suggest a spe- 
cific growth pattern as a function of cosmic time. According 
to these authors, BHs with M B h ~ 10 8 M G attained 50% of 
their mass by z = 2 and 90% of their mass by z = 1 (Mar- 
coni et al. 2004, Fig. 2). Smaller BHs grew slower at ear- 
lier times and many active BHs with Mbh ~ 10 7 M© are still 
growing today. Similar scenarios, under the general terminol- 
ogy of "anti-hierarchical growth of supermassive BHs", have 
been presented by Merloni (2004) and others. Those stud- 
ies assume that the growth rate of very massive BHs at high 
redshifts approached the Eddington limit. More recent stud- 
ies, e.g., by Volonteri et al. (2006) and Hopkins et al. (2006), 
focus on the importance of the Eddington ratio, L/Leaa, in de- 
termining BH evolution at all redshifts. This includes also hi- 
erarchical models for the evolution of the most massive BHs. 
All these models can be tested observationally by direct mea- 
surements of BH mass and accretion rate at high redshifts, 
provided high quality observations and reliable methods for 
determining Mbh are available. 

Current BH mass estimates are based on reverberation map- 
ping that provides a way to measure the emissivity-weighted 
size of the broad-line region (BLR) in type-I AGNs as a func- 
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tion of the optical continuum luminosity {XLx at 5 100A, here- 
after L5100; see Kaspi et al. 2000; Kaspi et al. 2005, hereafter 
K05 and references therein). This scaling has been used to 
obtain a "single-epoch" estimate of the BH mass by combin- 
ing the BLR size derived from L5100 with a measure of the 
gas velocity obtained from FWHM(Hj3) (hereafter "the H/3 
method"). There are obvious limitations to this method due 
to the somewhat vague definition of the BLR size, variable 
source luminosity, BLR geometry, and somewhat uncertain 
line widths. There are also questions regarding the exact slope 
of the BLR-size vs. the source luminosity (e.g. Bentz et al. 
2006). These translate to a factor of ~ 2 uncertainty on the 
derived masses. An additional uncertainty is associated with 
the limited luminosity range of the K05 sample and hence the 
need to extrapolate the relationship beyond its highest lumi- 
nosity end, at L$iqq~ 2 x 10 46 ergs s _1 (see, however, the new 
result of Kaspi et al. 2007). 

Other combinations of continuum luminosity and line 
widths have also been used, especially in the study of high- 
redshift sources (e.g., McLure & Dunlop 2004; Vestergaard 
2004; Vestergaard & Peterson 2006). These seem to be asso- 
ciated with a larger uncertainty on the derived masses espe- 
cially when the C IV X 1549 line width, in combination with 
XLx at ~ 1400 A (hereafter "the C IV X 1549 method") are 
used (e.g., Baskin and Laor 2005). A way to extend the 
Hj3 method to high redshifts is by high quality near-infrared 
spectroscopy. Measurements in a moderate size sample (29 
sources) of this type are presented in Shemmer et al. (2004, 
hereafter S04) who used H- and /T-band spectroscopy for ob- 
taining single-epoch BH masses and accretion rates in very 
high luminosity AGNs. More recent work of this type (a sam- 
ple of 9 sources with 1 .08 < z < 2.32) is reported in Sulentic 
et al. (2006). 

This paper presents BH masses and accretion rates obtained 
with the Hj3 method for a new sample of 15 sources at two 
redshift bands around z — 2.3 and z — 3.4. Together with our 
earlier work (S04) we can now use the Hj3 method to investi- 
gate the mass and accretion rate of high redshift sources over a 
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TABLE 1 



GNIRS Observation Log 


Object ID (SDSS J) 


2sys a 


ZSDSS b 


Band 


Obs. Date 


Exp. Time (hr) 


025438.37+002132.8 


2.456 


2.463 


H 


2005 Dec 7 


1.0 


083630.55+062044.8 


3.397 


3.397 


K 


2005 Nov 27 


1.0 


aac 1,11 nn i moicn c 

095141.33+013259.5 


2.411 


2.429 


H 


2006 Jan 24 


2.5 


100710 70+0421 19 1 


2.363 


2.363 


n 


2006 Mnr 26 


2.0 


101257.52+025933.2 


2.434 


2.434 


H 


2006 Feb 08 


0.5 


105511.99+020751.9 


3.391 


3.384 


K 


2005 Apr 21 


1.0 


113838.26-020607.2 


3.352 


3.343 


K 


2005 Apr 21 


1.0 


115111.20+034048.3 


2.337 


2.337 


H 


2006 Apr 16,25 


2.5 


115304.62+035951.5 


3.426 


3.432 


K 


2005 Apr 23 


0.7 


115935.64+042420.0 


3.451 


3.448 


K 


2005 Apr 23 


0.7 


125034.41-010510.5 


2.397 


2.397 


H 


2006 Apr 26,27 


2.7 


144245.66-024250.1 


2.356 


2.343 


H 


2006 Jul 19 


1.0 


153725.36-014650.3 


3.452 


3.452 


K 


2005 Apr 21 


0.7 


210258.21+002023.4 


3.328 


3.343 


K 


2006 Jul 19 


1.3 


210311.69-060059.4 


3.336 


3.336 


K 


2005 May 31 


0.7 



"Systemic redshift measured from the [O III] lines. 

b Redshift obtained from the SDSS archive, based on rest-frame UV emission lines. 



range of 2.5 dex in luminosity. In § 2 we describe the observa- 
tions and their analysis, and in § 3 we present the main results 
and a discussion of our findings including an assessment of 
the growth rate of high-redshift AGNs, and the evaluation of 
BH mass measurements obtained with the C IV X 1549 method 
as a replacement for the H/3 method at high redshifts. 

2. OBSERVATIONS, DATA REDUCTION, AND MASS 
DETERMINATION 

The sample described in this paper contains 15 high red- 
shift, high-luminosity AGNs that were selected by their lumi- 
nosity and redshift. The redshift range stems from the re- 
quirement to directly measure the H/3 line and the 5 100 A 
continuum from the ground in order to use the H/3 method, 
which we consider to provide the most reliable BH mass es- 
timates. The redshift ranges are dictated by the wavelengths 
of the H and K bands and are ^2.1-2.5 and ^3.2-3.5, re- 
spectively. The luminosity is dictated by the goal of going 
down the AGN luminosity function, starting from its top, and 
measuring masses and accretion rates of fainter and fainter 
sources. Given this, we chose sources that are 5-10 times less 
luminous than the mean luminosity in the S04 sample. For the 
basic sample we chose the Sloan Digital Sky Survey (SDSS; 
e.g., York et al. 2000) which contains flux calibrated data and 
hence luminosity estimates for all sources. 

Spectroscopic observations were obtained with the Gemini 
Near Infrared Spectrograph (GNIRS) on Gemini South under 
programs GS-2005B-Q-28, GS-2006A-Q-58 and GS-2005A- 
Q-5 1 . The long slit observations were acquired in the H or 
K bands depending on the redshift of the source. The slit 
width was 1" and the targets were nodded along the slit to 
obtain a good background subtraction. The 32 1/mm grating 
was used in all observations resulting in R ~ 640 and 850, 
AX ~ 1.46- 1.89 jum and AX ~ 1.83 -2.49 jlim for the H 
and K bands, respectively. Typical exposure times of the 
sub-integrations were 400 to 600 seconds. More details are 
given in Table 1 where we list all 15 sources and assign to 
each the systemic redshift measured from their [O III] A5007 
lines. These can be somewhat different from the SDSS red- 
shifts that were obtained by measuring rest-frame ultraviolet 
(UV) broad-emission lines and are also listed in the table. 

The reduction of the raw spectroscopic data was done using 
the GEMINI package in IRAF. The pipeline combines obser- 



vations from different nodded positions to obtain background- 
subtracted images, determines the wavelength calibration, 
registers the frames, and produces a final averaged image. 
Extraction of the spectra and flux calibration were performed 
using standard IRAF tasks. Special care was taken to correct 
for telluric absorption. This was done by observing early- 
type stars right before or after the science targets and at simi- 
lar air masses. 

To obtain a more accurate flux calibration of our spectra, H- 
and /f-band photometry was obtained for 13 of our sources 
using the ISPI detector on the Cerro Tololo Inter-American 
Observatory 4 m telescope on 2007 February 12. The data 
were reduced in the standard way using the XDIMSUN pack- 
age in IRAF 5 . The calibration was achieved using Two Mi- 
cron All Sky Survey (2MASS 6 ) stars in the 10' x 10' field of 
view of the instrument. Only stars with good quality flags 
in the 2MASS All-Sky Catalog of Point Sources (Cutri et al. 
2003) were used, resulting in a ~ 10% accuracy. All luminosi- 
ties listed below are based on these values except for the two 
sources that were not observed where we use the flux from 
the spectroscopy. The typical RMS difference between the 
two methods is about 15%. 

We applied a fitting procedure similar to the one described 
in Netzer & Trakhtenbrot (2007; hereafter NT07) to fit all 
spectra with the various components expected in this range. In 
short, we fit a linear continuum between rest-frame ~ 4700A 
and ~ 5100A and then a five-component Gaussian emission 
line model to the continuum-subtracted spectrum: two com- 
ponents for the broad H/3 line, one for the narrow H/3 line, 
and two for the [O III] AA4959, 5007 A lines. This fit serves to 
obtain a first estimate of FWHM(H/3). We then use the Boro- 
son & Green (1992) Fe II template, convolved with a single 
Gaussian with the above FWHM (i.e. the one obtained from 
the combination of the two broad components), to fit the Fe II 
emission complex over the range 4400^-650A. This fit, ex- 
tended to the entire wavelength range by using the template, is 
used to estimate the Fe II line contributions to the continuum 
bands and to improve the continuum definition. Using the ini- 

5 IRAF (Image Reduction and Analysis Facility) is distributed by the Na- 
tional Optical Astronomy Observatories, which are operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 

6 http ://w w w. ipac . caltech.edu/2mass 
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FIG. 1 . — GNIRS spectra of our sample of 15 high-redshift AGNs. The best-fit model (thick solid curve) in each panel is composed of a continuum component 
(dotted curve), an Fe II-emission complex (thin dot-dashed curve), and Hj3 and [O III] emission lines (thick dashed curve). 



tial Fe II template, we obtain a modified continuum-subtracted 
spectrum and repeat the line fitting process. A second itera- 
tion Fe II model is obtained and then subtracted from the spec- 
trum, and a final, five-component Gaussian fit of the Hj3 and 
[O III] lines is performed. In the final stage, the two BLR 
Gaussians are limited to l,500<FWHM<20,000km s _1 and 
the NLR components are forced to have the same FWHM for 
all three lines. This component is not allowed to exceed 1,200 
km s -1 , which is smaller than the FWHM([0 III]) values of 
some of the S04 sources (see Table 1 of Netzer et al. 2004) but 
appears to be adequate for the present sample. The monochro- 
matic luminosity at 5100A, L5100, is measured directly from 
the fitted continuum. The results of our line-fitting procedure 
are given in Fig. 1 where we show the observed spectrum, the 
model components, and the final fit for all 15 sources. 

The measurement of Mbh is a crucial point and requires 
some discussion. The K05 study suggests that the corre- 
lation of BLR size as deduced from the Hj3 lag with re- 
spect to the optical continuum is given by Rblr L a where 
a = 0.65 — 0.7. This slope is obtained from fitting 34 sources 
covering about four orders of magnitude in L5 ioo- More recent 
work by Bentz et al. (2006) used an improved stellar subtrac- 
tion procedure for several of the lower luminosity sources in 
the K05 sample which reduced L5100 f° r those cases. Using 26 
sources they obtained a = 0.52. A more recent study by the 
same group (Bentz et al. 2007) that includes four additional 
sources (but again not the entire K05 sample) with improved 
stellar subtraction gives a = 0.54 ± 0.04. 

Given the very large luminosity range (about a factor 
2 x 10 4 in L5100) °f the K05 sample, it is not at all clear that 
the slopes at low and high luminosities are the same. It is 
thus justified to use the higher luminosity sources to obtain 
the most appropriate slope for extrapolating to luminosities 
larger than those of the K05 sample. Similarly, we suggest 
to use only the lower luminosity sources in K05 (after includ- 
ing the Bentz et al. 2007 corrections) when looking for the 
best extrapolation to very low luminosities. This approach 
was used by NT07 in their study of the SDSS sample. Since 



most of the source analyzed here have extremely large L5100 



(10 



,45.2-47.3 



ergs s '), we chose to adopt the same approach 



and obtain a best slope by fitting only those sources in K05 
with L5ioo> 10 43 - 5 ergs s _1 . This gives the same expression 
as used by NT07, 



M BH = 1.05 x 10 8 



£5100 



10 46 ergss- 



0.65 



FWHM(HJ3) 
lO^ms- 1 



Mr. 



(1) 

For comparison, we also calculated M Bli and accretion rate 
using the expression given in Bentz et al. (2007 Eq. 2). This 
gives masses that are smaller by a factor of ~ 1 .2 for our low- 
est L5100 sources, and by a factor of ~ 1.9 for our highest 
luminosity AGNs. These alternative estimates are compared 
below with the results obtained by using Eq. 1 . 

The calculation of Lboi/^Edd (hereafter L/Lgdd) is based 
on a bolometric correction, fi, which is somewhat lu- 
minosity dependent, and for the sample in hand is of 
order 5-7 (NT07 and references therein). For most 
applications described here we assumed f L — 1 [i.e., 
V L Edd=7xL 510 o/(l-5 x 1O 38 M BH /M )]. Marconi et al. 
(2004) give a specific expression for fi that translates to bolo- 
metric correction factors of 5.4-6.4 for our sample (see also 
Richards et al. 2006 on this issue but note the different method 
of counting the IR luminosity in that paper). We have car- 
ried the analysis described below using those values but given 
the small range in fi, the conclusions hardly change. Ta- 
ble 2 shows the results obtained from the above measurements 
and fits by assuming a standard cosmology with Q.a = 0.7, 
Q. m = 0.3, and H = 70 km s _1 Mpc" 1 . 

All the 29 spectra of the S04 sample were refitted using 
a similar procedure and BH masses were recalculated us- 
ing Eq. 1 (which is slightly different from the one used in 
S04). Some of the S04 data are of poorer quality compared 
to the new GNIRS observations which resulted in larger un- 
certainties. In addition, slight mis-centering of the Hj3 line 
in the H or K bands forced us, in several cases, to perform 
the Fe II fit on the Fe II complex longword of the Hj3 line. 
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TABLE 2 

Observed and Derived Properties 



Object ID (SDSS J) z sys logLjlOO log AL A (1450A) FWHM(Hj8) FWHM(C IV) log M BH logL/L E dd %row a / 

(ergss" 1 ) (ergss" 1 ) (kms" 1 ) (kms" 1 ) (M ) funiverse 



025438.37+002132.8 


2.456 


45.85 


45.93 


083630.55+062044.8 


3.397 


45.53 


46.12 


095141.334 


-013259.5 


2.411 


45.55 


45.90 


100710.704 


-042119.1 


2.363 


45.17 


45.67 


101257.524 


-025933.2 


2.434 


45.73 


45.95 


105511.994 


-020751.9 


3.391 


45.70 


46.24 


113838.26- 


-020607.2 


3.352 


45.79 


46.18 


115111.204 


-034048.3 


2.337 


45.58 


45.58 


115304.624 


-035951.5 


3.426 


46.04 


46.37 


115935.644 


-042420.0 


3.451 


45.92 


46.43 


125034.41- 


-010510.5 


2.397 


45.41 


45.71 


144245.66- 


-024250.1 


2.356 


46.03 


45.90 


153725.36- 


-014650.3 


3.452 


45.98 


46.44 


210258.214 


-002023.4 


3.328 


45.79 


45.86 


210311.69- 


-060059.4 


3.336 


46.30 


46.24 


"Assumin; 


gJJ=0.1,/ L 


= 7, -M S eec 


l = 1O 4 M , and / ac ti v e 


= 1. 



Another difference from the S04 procedure is the inclusion 
of the narrow Hj3 component in the fits. As a result, sev- 
eral of the newly measured FWHM(H/3) values are some- 
what larger than those found by S04. Given all this, we 
consider the newly fitted FWHM(Hj3) values to be more re- 
liable than those presented in S04, although the differences 
are small and the main change is the inclusion of an es- 
timate on the uncertainty of the FWHM(Hj3) measurement 
(see below). Notable exceptions are six cases (UM632, 
2QZJ231456.8-280102, [HB89] 2254+024, SBS 1425+606, 
UM642 and SDSS J024933.42-083454.4) where the differ- 
ence in Mbh is of order ^2-3 due to the change in the mea- 
sured FWHM(Hj3). The total sample consists of 44 sources 
covering the luminosity range 10 45 2 < XL X (5100A) < 10 47 3 
ergs s _1 . This represents the AGNs luminosity function from 
about a factor of 2 below the top to a factor of ~ 200 below 
the top. 

The main uncertainties on the measured masses are due to 
uncertainties in FWHM(Hj3). To estimate this, we divided 
all fits into three categories reflecting their quality. Most of 
the sources show symmetrical lines and the Fe II complex is 
easy to model and deconvolve. The assigned uncertainty on 
FWHM(Hj3) in this case is 10%. Some sources have adequate 
S/N ratios yet the broad line profiles are somewhat irregular 
and the FWHM is more difficult to constrain. The uncer- 
tainty in this case is estimated to be 20%. Finally, in those 
cases showing asymmetric profiles, difficult to model Fe II 
lines and poorer S/N, we assigned an uncertainty of 30% on 
FWHM(Hj3). While This procedure is somewhat subjective, 
we have no better way to quantify the fitting process. We con- 
sider those uncertainties conservative and note that they trans- 
late to relative errors on the mass determination of 20-60%. 
None of the new GNIRS objects is assigned the highest uncer- 
tainty of 30% but six of the S04 sources fall into this category. 

We also used the SDSS spectra to measure AL A (1450A) and 
FWHM(C IV A 1549) for the 15 new sources. The luminosity 
is available directly from the observed continuum flux and 
the FWHM is obtained by fitting two Gaussians to the line 
profile. Similar information is available in S04 for 27 of the 
29 sources in their sample. 

To summarize, our sample contains almost all sources at 
z > 2 where Mbh was obtained using the most robust and reli- 
able method (the H/3 method). Much larger uncertainties are 



4164 


4753 


9.162 


-0.64+0.08 


0.8 


3950 


5878 


8.909 


-0.71+0.15 


1.3 


4297 


5289 


8.997 


-0.78+0.15 


1.1 


5516 


5495 


8.96 


-1.13+0.15 


2.4 


3892 


5862 


9.029 


-0.62+0.08 


0.8 


5424 


5476 


9.294 


-0.93+0.08 


2.4 


4562 


6098 


9.271 


-0.74+0.15 


1.5 


5146 


2860 


9.171 


-0.92+0.08 


1.5 


5521 


1773 


9.529 


-0.82+0.08 


2.0 


5557 


4160 


9.460 


-0.89+0.15 


2.2 


5149 


4234 


9.061 


-0.98+0.08 


1.8 


3661 


3277 


9.166 


-0.47+0.08 


0.6 


3656 


5650 


9.133 


-0.49+0.08 


0.9 


7198 


2355 


9.599 


-1.14+0.08 


4.0 


6075 


4951 


9.785 


-0.81+0.08 


2.0 



associated with the C IV A 1549 method (see detailed discus- 
sion below), the one used in almost all other z > 2 studies. 
Thus, we believe that our data set is the most suitable and 
most accurate to address the issues of BH growth and the dis- 
tribution of Mbh and L/L^dd at those redshifts. 

3. DISCUSSION 
3.1. Luminosity and L /LEdd in High-Redshift AGNs 

We have looked for correlations between L/L E ^d and var- 
ious other properties of the 44 high-redshift sources in our 
sample. The correlation with L5100 is shown in Fig. 2a. The 
two sub-samples are shown with different symbols emphasiz- 
ing the difference in the mean luminosity. Standard regression 
analysis indicates a significant correlation with a Pearson cor- 
relation coefficient is 0.45 (p — 2 x 10~ 3 ) and a Spearman- 
rank correlation coefficient of 0.46 (p = 2 x 10~ 3 ). We also 
checked the correlation of L/L^ with Mbh- The diagram is 
shown in Fig. 2b and exhibits a large scatter and no apparent 
correlation. In particular, the scatter in L/LEdd (about a factor 
of 10) is similar for BHs of all masses. 

We have repeated the analysis using, this time, the Bentz et 
al. (2007) expression for estimating Mbh- The correlation of 
L/L E< id with L5100 is even stronger and, again, there is no cor- 
relation of L/L E dd with M B h- As explained, all M B h values 
calculated in this way are smaller by factors of 1.2-1.9, de- 
pending on their luminosity. As a result, all values of L/L^ 
are larger by similar factors. This results in some sources with 
L/LEdd^ 3. We suspect that the extremely large accretion 
rates may not be physical but given the method uncertainty, 
and the extrapolation beyond the K05 luminosity range, we 
cannot rule them out. On the other hand, the deduced Mbh in 
this case is smaller, on average by a factor 1.5, which may be 
more consistent with the lack of very high mass BHs at low 
redshifts (but note that even the Bentz et al. relationship gives 
several cases with Mbh > 1O 1O M ). 

The slopes of the above correlations depend on the statis- 
tical method used and are not too different from the one ex- 
pected from a case where FWHM(Hj3) is independent of L5 100 
(Eq. 1). Indeed, there is no correlation between FWHM(Hj3) 
and L5100 in our sample. This point requires some explana- 
tion. Single-epoch mass determination provides reliable BH 
mass estimates only because of a (yet to be explained) scal- 
ing of the BLR size with source luminosity (discovered by 
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FIG. 2. — i/iEdd vs - ^5100 andMBH (fight) for the entire sample of 44 
luminous, high-redshift AGNs. GNIRS (S04) sources are marked with open 
(full) symbols. The dashed and dotted lines mark two best-fit lines obtained 
with two, somewhat different statistical methods with slopes 0.29 and 0.43. 



reverberation mapping) and the virial motion of the BLR gas. 
Given this scaling, a complete and unbiased sample must also 
show some dependence of the mean gas velocity on source 
luminosity, depending on the distribution of Mbh in the sam- 
ple. This is not observed in the sample at hand. It may re- 
flect the incompleteness of the sample, its small size, or the 
real M Bli distribution. Thus, the above correlations are not, 
by themselves, very important. The more significant finding 
is the presence of a large number of very massive BHs, at 
high redshifts, with L/L^dd considerably smaller than unity. 
In this respect, there is no difference between the two redshift 
groups presented here. As discussed below, this is relevant to 
the question of BH growth in the early universe. 

The recent work of Kollmeier et al. (2006; hereafter 
K06) includes a systematic study of BH masses and accre- 
tion rates in a sample of 407 AGNs covering the redshift 
range of 0.3-4. These authors suggest a very narrow range 
of L/Ledd, at all luminosities and redshifts, consistent with 
log (L/Ledd)=— 0.6 ± 0.3. The range is even smaller (0.28 
dex) for 131 high luminosity, high redshift (z > 1.2 by their 
definition) sources that are more relevant to the present dis- 
cussion. The paper suggests that the intrinsic distribution in 
L/Lgdd is even narrower (0.24 dex for the high-z high-L sub- 
group) and much of the observed scatter is due to uncertainties 
in BH mass determination and bolometric correction. 

Our sample contains a similar number of sources to K06 
in the 2.1-3.5 redshift range. Moreover, the number of K06 
sources in the redshift and Mbh range of our sample is much 
smaller (only 19 compared with our 44). Thus our sample is 
more suitable, in terms of number of sources, to address the 
issue of the L/L E dd distribution in the population of high red- 
shift large BH mass AGNs. We find a broader range in L/LEdd 
compared with the various sub-groups presented in K06. For 
example, 90% of the 44 sources are found in the accretion rate 
interval 0.08 <L/Lgdd< 1 -5. It is not entirely clear what is the 
source of the difference between our results and those of K06. 
It may be related to the large uncertainty, and perhaps even a 
bias, in the method they used to determine Mbh (most masses 
in their high-z high-L group and all masses for z > 2 sources 
were determined with the C IV X 1549 method; see also com- 
ment on L/Lgdd measured this way in § 3.4). The K06 sample 
is flux limited and thus more complete than ours. However, 
as explained, the number of sources in individual mass and 
redshift bins are extremely small. These conclusions do not 
change when using the Bentz et al. (2007) M B h estimate since 



the discrepancy is mostly due to the range in L/L^dd which is 
even larger when this method is used. 

To illustrate the above points, we show in Fig. 3a one of 
the K06 histograms (their Fig. 11, panel with 2 < z < 3 and 
Mbh= 10 9 ~ 10 ) alongside the 34 objects we observed within 
a similar mass and redshift range. A visual inspection shows 
the broader distribution of L/L E dd m our sample (note that the 
K06 distribution should be shifted to the left by about 0. 1 dex 
due to the different bolometric correction used in their work). 

A more quantitative test can be made by comparing our 
sample with the entire high-luminosity high-redshift sub- 
sample of 131 objects from K06 (Table 1), despite of the dif- 
ferent range in z and Mbh- The overlap between the two is 
large but the mean luminosity and BH mass in our sample 
is somewhat larger. According to K06, log(L/L E dd) for this 
group is well fitted with a log-normal distribution which is 
centered at —0.52 and has a measured (model) standard devi- 
ation of 0.28 (0.24). As explained, the mean should be shifted 
to —0.62 to allow for the somewhat larger bolometric correc- 
tion used by K06. Our sample of 44 sources shows a similar 
mean (—0.56) but a considerably larger scatter of <~ 0.35. It 
also shows a positive skewness of <~ 0.47, in comparison with 
the K06 value of —0.02. This implies that sources with larger 
values of L/L E dd are more abundant in our sample. All these 
points are illustrated in Fig. 3b where we compare the his- 
togram of our measured values of L/L E dd with the favored 
K06 distribution for the 131 high-z high-L sources. 

We attempted to verify the null hypothesis that the val- 
ues of L/Ledd in our sample are indeed drawn from a log- 
normal distribution. A two-tailed Kolmogorov-Smirnov test 
over a large range of possible mean and standard deviation 
gives inconclusive results. The highest probability case has 
a mean log(L/L E dd)= —0.59 and a = 0.37. The /rvalue of 
this case is 0.92 (i.e. the probability of such a log-normal 
distribution in our sample is 92%). We also find that a K06 
distribution (log(L/L E dd)= —0.62 ±0.24 given our bolomet- 
ric correction) can be rejected at the 92% level. Finally, we 
tested the suggestion of a very large mean L/L E dd, close to 
unity, in a sample like ours. A series of Kolmogorov-Smirnov 
tests show that a log-normal distribution of any width can be 
rejected at the 99% confidence level for all cases where the 
mean log(L/L E dd) is —0.25 or larger. As shown in § 3.2, this 
is relevant for the comparison with BH evolution models. 

To conclude, our sample of high luminosity high redshift 
sources seems to be characterized by a similar mean but a 
somewhat broader distribution in L/L Ed d compared with the 
similar properties sub-samples in K06. Clearly some of the 
differences may be related to the way we chose our sample 
and in particular the fact that it is not a real flux limited sam- 
ple. This prevents us from reaching firm conclusions about 
the source of the difference at this stage. However, the method 
used here is, in our opinion, the best for measuring Mbh an d 
is preferable to the K06 method (see § 3.4). Given all this, 
we suggest that the distributions of our 44 sources shown in 
Fig. 3 represent the AGN population, over this redshift and 
M B h range, in the best way. The large range in L/L E dd (a fac- 
tor of ~ 20) is probably real and is likely to represent correctly 
the intrinsic properties of the AGN population analyzed here. 

3.2. Comparison with BH Evolution Models 

The cosmic evolution of BHs can be modeled by using 
large, recently observed AGN samples. This requires a com- 
bination of the X-ray and optical luminosity functions (LFs) 
since the first is more complete yet the latter probes much 
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FIG. 3. — Left: Histograms of L/Le<m values as obtained from K06 (Fig. 11, 
for 2 < z < 3 sources with Mbh = 1O 9 ~ 1O M , dotted line) and the present 
work (solid line, Mbh = 10 8 - 9 ~ 10 M Q , using f L = 7). Note that we directly 
copied the data from the K06 paper, thus, the somewhat different bolometric 
correction used by those authors would shift their histogram to the left by 
about 0.1 dex. Right: A comparison of the log(L/Lrjdd) distribution in our 
sample with the favored K06 model (dashed line). 



deeper at high redshifts. Of the numerous papers discussing 
such models we focus on the works of Marconi et al. (2004), 
Merloni (2004), Hopkins et al. (2006) and Volonteri et al. 
(2006). 

Marconi et al. (2004) and Merloni (2004) address the 
growth rate of active BHs assuming a constant accretion rate, 
i/Z/Edd=l- Both papers suggest faster growth rate of very 
massive BHs at earlier epochs in order to fit the observed LFs. 
The more recent studies of Hopkins et al. (2006) and Volon- 
teri et al. (2006) include estimates of the typical L/L^dd re- 
quired to match the redshift-dependent LFs. According to the 
modeling of Hopkins et al., history of accretion and the ob- 
served properties are not identical since the low accretion rate 
phase is hardly observed due to obscuration. In particular, the 
large mass very high luminosity AGNs are observable only 
close to their peak luminosity and accretion rate. This results 
in a narrow accretion rate range around L/L^a=\, consistent 
with the Marconi et al. (2004) assumption. Volonteri et al. 
(2006) used a combination of the Hopkins et al. results with 
theoretical merger rates expected in cold dark matter scenar- 
ios. The accretion rate in their models is changing with red- 
shift with 0.3 <L/L E dd< 1 for 3 < z < 6. 

The typical uncertainty on the values of L/L E dd m our sam- 
ple is a factor of ~ 2 and hence the mean value is not sig- 
nificantly different from the above predictions with their own 
large uncertainties. However, the trend seems to be different 
and points in a different direction. In particular, the theoretical 
requirement of L/LEdd— 1 at high redshift and large BH mass 
is not seen in many of our sources. In fact, the range of L/LEdd 
covered by our very high luminosity AGNs, from about 0.07 
to about 1.7 (with the Bentz et al. 2007 expression the range 
is 0.08-3.2), is not very different from the range typically ob- 
served in low-redshift AGN samples that are 2-3 orders of 
magnitude less luminous. Thus, it seems that some large mass 
high-redshift active BHs are observed far from their peak lu- 
minosity phase which presents a challenge to current theo- 
retical models. Unfortunately, our sample does not go deep 
enough to search for even lower L/L E dd sources at high red- 



shifts and the K06 data does not help either since all their z > 2 
mass measurements are based on the C IV A 1549 method. 

Finally, we mention several recent publications that attempt 
to measure mass and accretion rates in even higher redshift 
sources. Jiang et al. (2006) presented Spitzer photometry 
of 13 AGNs at z ~ 6. The BH mass and accretion rate of 
four of the sources were estimated by the C IV X 1549 method, 
and the M Bii and L/LEdd values are given in their Table 3. 
All four BH masses are very large, ~ 5 x 10 9 M and the 
L/Lem range is 0.5-1. However, the bolometric correction 
used by these authors is based on integrating over the en- 
tire spectrum, from hard-X-rays to the far-infrared. This in- 
volves double-counting since almost half of this emission is 
due to reprocessing of the primary continuum (accretion disk 
and X-ray source) by the dusty medium near the BH. Thus, 
the "true" L/L^dd values are smaller than those listed in their 
paper by a factor of ^2. Moreover, the BH mass estimates 
are rather uncertain because of the use of the C IV X 1549 
method (see § 3.4). Results of similar quality on a handful 
of additional AGNs at z ~ 6 have also appeared recently (e.g., 
Willott et al. 2003; Jiang et al. 2007; Kurk et al. 2007). 

3.3. The Growth Rate and Growth Time of High Mass, 
High-Redshift BHs 

The combination of the newly measured BH masses and 
accretion rates can be used to estimate the growth times of 
massive BHs at z ~ 2.3 and z ~ 3.4. We follow the procedure 
outlined in NT07 (their Eq. 6) to calculate f grow assuming 



*grow 



?Edd 



77/(1-77) 



■log 



(M BH \ 



1 



(2) 



/i^5100/i-Edd \M see dJ /active 

where ? E dd = 3.8 x 10 8 yr for cosmic abundance, 77 is the ac- 
cretion efficiency, and /active is the duty cycle (the fractional 
activity time) of the BH. The growth time, f grow , is most sen- 
sitive to 77 and /active since both can vary by large factors. 
The seed BH mass, M see d, may also change over a large range 
but /g ro w is not very sensitive to this change. A small seed 
BH (10 2 - 10 3 M©) can be the result of population III stars, 
at z ~ 20, while larger seed BHs (10 4 - 1O 6 M ) can be due 
to direct collapse at lower redshifts (Begelman, Volonteri, & 
Rees, 2006). Finally, f grow is less sensitive to f L because of the 
limited expected range of this factor in the sample at hand. 

As discussed in several publications, most recently by King 
& Pringle (2006), the value of 77 depends on the direction 
and magnitude of the BH spin and the angular momentum of 
the accretion disk. The highest value is ~ 0.4 and the low- 
est values (<~ 0.04) is obtained for retrograde accretion with 
a BH spin parameter of a = — 1. More typical values, that 
reflect several spin-up and spin-down episodes with small- 
mass disks (see King & Pringle 2006), give 77 in the range 
0.05-0.08. The values of f grow calculated here (Table 2) are for 
the case of 77 = 0.1, M sced = 10 4 M Q , / ac tive = 1, and f L = 7. 
We list them as fg row / Universe where /universe is calculated at 
the given redshift of the source (using our adopted cosmology 
from § 2). 

The computed values of f grow exceed the age of the uni- 
verse in most of the lower L/Leaa sources. While somewhat 
smaller values of f gro wA universe are indeed likely, given the 
expected range in 77, we consider the numbers given in Ta- 
ble 2 as conservative lower limits since / ac tive is likely to be 
much smaller than unity. For example, general galaxy and 
BH growth considerations predict / ac tive ~ 0.1. This can be 
estimated from number counts of AGNs and galaxies at large 
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redshifts, from estimates of typical time scales for powerful 
starbursts, and from models of BH and galaxy evolution. An- 
other estimate is obtained by assuming that the ratio of BH 
mass to the host galaxy mass is similar to the typical value 
observed in the local universe (about 1/700). This translates, 
for most objects in our sample, to a host mass of 1O 12 M or 
larger. To build up such a large stellar mass would require 
continuous star formation, from very early times, at a rate of 
~ 300M Q /yr, or a much higher star formation rate with a star- 
formation duty cycle which is less than unity. Given the much 
faster growth of BHs by accretion, their duty cycles would 
be considerably shorter. While some very different scenarios 
cannot be excluded on the basis of current observations, they 
are definitely not in accord with present day galaxy-formation 
models. Anti-hierarchal BH growth models would give faster 
BH growth at high redshifts with a larger BH-to-galaxy mass 
ratio. It remains to be seen whether any such model is consis- 
tent with /active ~ 1 all the way to redshift 3 for BHs like the 
ones in our sample. 

Given our more plausible assumption of /active -C 1 for 
z = 2 — 3 very large mass BHs, all 15 new sources presented 
here, and several of the remaining more luminous objects, 
require too much time to grow to their observed size. This 
problem was not noted in earlier works since most of them 
assumed L/L^^ 1 and hence f gr0 w an order of magnitude 
shorter. The problem is most severe for the lowest accretion 
rate sources. The higher accretion rate sources (most of the 
29 objects from S04) have just enough time to explain their 
mass if / ac tive ~ 1 • Thus, our new measurements, and general 
considerations, indicate that some very massive BHs had just 
enough time to grow to their measured size at redshifts <~ 2.3 
and ~ 3.4. However, for all sources with L/Ledd ~ 0.3, there 
was not enough time at redshift ~ 3.4, and even at redshift 
~ 2.3, to grow to their observed size. This problem is more 
severe if we adopt the Hopkins et al. (2006) set of models 
where much of the BH growth is during times where L/^Edd 
is smaller than the time where the object is not obscured. 

We conclude that a significant fraction of the sources in our 
sample must have had at least one previous episode of faster 
growth, probably with L/Ledd^ 1 and at z £ 2 — 3, in order to 
explain their BH mass. Given this, many extremely large mass 
BHs at redshifts 2 and 3 are in the process of their second and 
perhaps even third or fourth episode of activity. This does not 
necessarily apply for small-mass BHs at those redshifts that 
may have a different growth pattern. 

3.4. C IV A 1549 as a BH-Mass Indicator 

As noted above, the C IV A 1549 method for estimating BH 
masses is problematic since C IV X 1549 lines are known to 
show unusual profiles in many cases, including an enhanced 
blue wing and a large wavelength shift relative to the sys- 
temic velocity. This issue has been discussed, extensively, 
in numerous papers. The work of Baskin & Laor (2005) 
clearly shows the differences between the H/3 and C IV A 1549 
line profiles and the problematics of FWHM(C IV A 1549) as a 
virial-motion indicator (see also the recent Shang et al. 2007 
work). However, according to Vestergaard & Peterson (2006), 
a proper calibration of FWHM(C IV A 1549) and the UV con- 
tinuum provides an adequate replacement for the Hj3 method. 
K06 adopted this view and used the C IV A 1549 method to 
estimate BH masses in many high-redshift sources. 

Our combined sample includes 44 extremely luminous 
AGNs with FWHM(H/3) and L 5100 values for all. We have 
compiled AL A (1450A) and FWHM(Civ A 1549) for 42 of 
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FIG. 4. — BH mass calculated with the Hj3 method vs. the one obtained 
from the combination of FWHM(C IV A 1549) and the UV continuum. Sym- 
bols are as in Fig. 2. The dashed line marks a 1 : 1 correspondence, to guide 
the eye. 



these and are thus in a position to compare the H/3 and 
Civ A 1549 methods in the high luminosity range. We first 
checked L5100 vs. AL^(1450A) for all sources where we 
have independent luminosities for both (i.e., excluding the six 
sources from S04 for which optical luminosities were esti- 
mated from their UV luminosities). We find a very significant, 
low scatter correlation which is consistent with L v oc v~ 5 . 
The diagram (not shown here) suggests that the two estimates 
of the source luminosity can be used, interchangeably, with 
a different normalization, for estimating the BLR size. We 
then obtained BH masses by using the two methods; the H/3 
method (Eq. 1) and the C IV A 1549 method with the expres- 
sion given in Vestergaard & Peterson (2006). The two mass 
estimates are shown in Fig. 4. The figure is a complete scatter 
diagram with no apparent correlation or trend. While the me- 
dian (0.89) and the mean (1.49) M BH (CrvAl549)/M BH (H/3) 
ratio are not far from unity, the scatter is about ±0.3 dex. 

Repeating the same analysis by using the Bentz et al. 
(2007) expression, gives also a scatter diagram with no sig- 
nificant correlation. In this one, the median and the mean for 
M BH (C IV A 1549)/M BH (H/3) are 1 .3 and 2.2, respectively and 
the scatter is about 0.4 dex. This may provide a slight indica- 
tion that the slope of a = 0.65 used in our work is preferable 
to the one used by Bentz et al. (2007) since, using this slope, 
the agreement between the median M B h(C IV A 1549) and the 
median M B h(H/3) is improved. 

The lack of correlation between M BH (H/3) and 
M bh (CivA1549) suggests that, while BH mass esti- 
mates based on the C IV A 1549 line properties roughly follow 
the increase in source luminosity, their use is rather limited 
and highly uncertain, in the high luminosity and/or high 
redshift range. This conclusion is independent of the slope 
used to obtain Af BH (H/3). As argued earlier, some of the 
differences found between our L/Leaa distribution and the 
one presented in K06 are likely to be the result of their use 
of the Civ A 1549 method. In fact, we notice a complete 
lack of correlation between the L/^Edd values derived with 
the H/3 method and those derived with the C IV A 1549 
method. Given all this we suspect that the unexpected trend 
found by K06 when comparing the various line-based mass 
measurements (see their Fig. 4 and related discussion) is due 
to the problematic C IV A 1549 method. 

4. CONCLUSIONS 
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We have presented new, high quality observations of the Hj3 
region in 15 high luminosity, high-redshift AGNs. We used 
the new data, in combination with our earlier observations of 
29 sources, to compile a sample of 44 high luminosity AGNs 
where the BH mass estimates are all based on the H/3 method. 
The main findings are: 

1 . Our sample is the largest of its type, yet, and can there- 
fore be used to investigate, in a more reliable way, the 
important correlations between source luminosity, BH 
mass and accretion rate. Such measurements at high 
redshifts can be combined with similar measurements 
at z < 0.75 to follow BH growth and AGN evolution 
through time. 

2. There is a significant correlation between L5100 and 
i/^Edd but no correlation between M B h and L/Ledd- 
These results do not depend on the exact slope used to 
derive Mbh- 

3. Assuming the distribution in L/Leu found here repre- 
sents the high-redshift high luminosity AGN popula- 
tion, we find a significant fraction of sources with small 
(< 0.2) L/Lgdd- Current theoretical models predict only 
very few such sources among high-redshift AGNs. 

4. Low accretion rate BHs represent a real challenge to 
BH growth scenarios. A possible way out is to as- 
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